Maximal Random Packing of Spheres in Menger Sponge 
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Abstract 

Random packing of spheres inside fractal collectors of dimension 2 < d < 3 is studied numerically 
using Random Sequential Adsorption algorithm. The paper focuses mainly on the measurement 
of maximal random packing ratio. Additionally, scaling properties of density autocorrelations in 
the obtained packing are analyzed. The RSA kinetics coefficients are also measured. Obtained 
results allow to test phenomenological relation between maximal random coverage ratio and col- 
lector dimension. Additionally, performed simulations together with previously obtained results 
confirm that, in general, the known dimensional relations are obeyed by systems having non-integer 
dimension at least for d < 3. 



* michal.ciesla@uj.cdu.pl 
t |ncbarbas@cyf- kr.edu.pl 



1 



I. INTRODUCTION 



Hard spheres packing has a reach research history, especially in maths but also in natural 
science, due to it's potential application. Attention has been mostly paid to maximal packing 
as it has many important implications, for example in crystallography The fundamental 
property defining packing is its density jl, 2|: 



Omax{d) = P'V(d), 

where d is dimension of space, p is a number of spheres in a unit volume, and 



(1) 



v(d) 



r(l + d/2) 



(2) 



is a volume of d dimensional sphere having radius r . is the Euler gamma function. 

The problem of maximal spheres packing is trivial in one dimension, simple in two di- 
mensions and only recently has been proved in three dimensional space [3|. For higher 
dimensions, it still remains open. 

Presented work, however, focuses on maximal random packing. It is derived from col- 
loids study, namely, modelling of irreversible adsorption processes, in which molecules are 
randomly placed on a collector surface. Once adsorbed, a molecule retains its position on a 
surface. Such adsorption layer grows till there is available uncovered place large enough to 
place subsequent molecules. As the simplest and most commonly used model of adsorbed 
molecule is sphere, the properties of adsorption layers are directly connected with random 
maximal packing of spheres. Here, the only analytically solved case is a one- dimensional 
problem, also known also as car parking problem, for which 9 max (l) = 0.748... j^]. Solutions 
in higher dimensions were found only numerically, however, it is worth to notice that there is 
a number of analytically solvable models for lattice-like collectors e.g. 5|, [6|. For flat, homo- 
geneous and continuous collectors, the most extensively studied case is a two-dimensional 



problem due to its potential application in chemistry and material science [7-lll|. The 
analysis of hard spheres random packing for 2 < d < 6, nicely reviewed in [lO], provides 
phenomenological relation, further slightly corrected later in 12fl , between maximal random 
coverage ratio and collector dimension: 

1 + C\d + c 2 d 2 



"max (d) 



2 d 



(3) 



with c\ = 0.416 and C2 = 0.086. There are also some characteristics, which do not depend 
qualitatively on dimension. For example, the two-point density correlation function for 
maximal random coverages is known to have super-exponential decay with growing distance 



13| 



C(r) 



(r/2r )-l 



T(r) \\n(r/(2r ) - 1 
and logarithmic singularity when particles are in touch 



for r — > 00, 



(4) 
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C(r) ~ In 



2r n 



1 ) for r ->■ (2r ) + and 9 -> 9 V 



(5) 



Aim of present work is to check validity of the above relations for adsorption on fractal 
collectors with structure similar to the Menger sponge. We decided to focus on this kind 
of collectors because its structure remains zeolites, which were extensively studied in the 
context of adsorption, e.g.jjjj]. 

Adsorption experiments using fractal-like collectors were made by Kinge et al. 



151 ] . though 



preparation of such surfaces is quite complicated. On the other hand, t 



rere are number of 



natural porous media when adsorption plays an important role there [16l. Il71j . For example, 
fractal-like structure of corals helps them to catch plankton effectively jis|. Adsorption on 
fractal collectors might also be applied in environmental protection in desi gnin g effective 



filters 



19[. Other studies on fractal collectors investigate diffusion properties 20|, |2l| or ad 



sorption on rough surfaces |22l. |23|. On the other hand, fractal structures were also observed 
as a result of adsorption [24|. However, the only similar study concerning adsorption on 
fractal collectors is limited to collectors of D < 2 12 1. 

The most common method of obtaining numerically the maximal random coverages uses 
Random Sequential Adsorption (RSA) algorithm [?| described in detail in SecJTTJ The ki- 
netics of RSA simulation has also been of our interest because it affects 9 max (d) estimation. 
Moreover, it was observed that for sufficiently long simulation time, the coverage ratio 9{t) 
scales with collector's dimension {7, 8]: 



9 max — 9(i) ~ t 



-l/d 



(6) 



Here 9 r . 



9{t — > 00). The above relation is all the more important as it is also valid for 



other molecules, more complex than simple spheres 
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II. MODEL 



Maximal random coverages are generated using RSA algorithm, which has been success- 
fully applied to study colloids. It is based on independent, repeated attempts of adding 
sphere to a covering film. The numerical procedure is carried out in the following steps: 

i: a virtual sphere is created with its centre position on a collector chosen randomly according 
to the uniform probability distribution; 

ii: an overlapping test is performed for previously adsorbed nearest neighbors of the virtual 
particle. The test checks if surface-to-surface distance between spheres is greater than 
zero; 

iii: if there is no overlap the virtual particle is irreversibly adsorbed and added to an existing 
covering layer. Its position does not change during further calculations; 

iiii: if there is an overlap the virtual sphere is removed and abandoned. 

Attempts are repeated iteratively. Their number is typically expressed in a dimensionless 
time units: 

to = N^, (7) 

where N is a number of attempts and Sb denotes a volume of a single sphere projection 
on a collector (note that only the center of a sphere has to be on collector) and Sq is 
a collector volume, assuming its integer dimension. Here, in D < 3 the volume of any 
collector is formally zero. On the other hand, any given value of Sc is only a scale factor 
of t . Therefore, for the purposes of our study, we assumed Sq = (100 ■ r ) 3 , which is the 
volume of an overlapping cube. Collectors were modelled as a successive iteration of a given 
fractal. They were constructed starting with a single cube, then dividing it into 27 subcubes 
and finally removing some of them. Then the procedure was repeated iteratively for each of 
the remaining subcubes. The real fractal is achieved after an infinite number of iterations. 
Its dimension is given by relation 

where p is a number of remaining subcubes in each iteration and s is a change of scale. In 
our case, s = 3. For example for Menger sponge p = 20 so D = ln(20)/ln(3) ~ 2.727. 



With each of the 27 subcubes marked with its Cartesian coordinates, the fractal can be fully 
described by coordinates of removed or retained subcubes. For convenience, basic properties 
of fractal collectors used in this study were collected in Tab]!] For each collector type 100 



removed subcubes l s< iteration of 

fractal dimension d 
coordinates a collector 



all but 113, 131, 222, 
311, 333 

all but 111, 113, 131, 
133, 222, 311, 313, 331, 
333 



1.4650 



2.0 



X 



112, 121, 123, 132, 211, 
213, 222, 231, 233, 311, 



313, 331, 333 



2.4022 




122, 211, 212, 213, 221, 



223, 231, 232, 233, 322 



2.5789 




122, 212, 221, 222, 223, 
232, 322 



2.7268 




211, 213, 231, 233 



2.8540 




222 



2.9656 




TABLE I: Definitions and fractal dimensions of collectors used for adsorption modelling. 
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independent numerical experiments have been performed. Each of them last 10 4 to- 

Typically coverage ratio is estimated by a number of adsorbed molecules. However, this 
approximation can be misleading when collector dimension differs from adsorbate dimension. 
Therefore, by analogy with integer dimension cases, covered area is limited to a cross-section 
of collector and adsorbate sphere. For this reason, the coverage ratio of a given adsorption 
layer is determined using random sampling of collector points and checking whether they 
are covered by any sphere adsorbed earlier: 

n c (t) 



6{t) 



n 



(9) 



where n c (t) is a number of covered points after a dimensionless time t; n denotes total 
number of sampled points. Here n = 10 6 , which provides statistical error at the level of 
0.1%. 

III. RESULTS AND DISCUSSION 

In performed simulations a real fractal was approximated by a its finite iteration. To 
check the level of possible systematic error introduced by such approximation we plotted 
maximal random coverage as a function of fractal iteration. Results presented in Fig{TJ show 
that measured values stabilize around 7 th iteration. Therefore, for the purpose of this study 
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FIG. 1: The dependence of maximal random coverage ratio on fractal iteration for all 

studied fractals. 



,we assumed that results obtained for 9 iteration math the true ones. 



A. Kinetics of the RSA 



The maximal random coverage is achieved after an infinite simulation time. Therefore to 
measure its value in finite simulations, the RSA kinetics model has to be used. The analysis 
of earlier numerical experiments shows that RSA for spheres in integer and fractal dimensions 



obeys the Feder's law ([6]) jlOl . l^ . The most convenient way to determine exponent in ([6]) 



is to find the time derivative of coverage ratio: 

fw~r- (10) 

where a = —1/d. Results presented in Fig J2] show that Feder's law slightly underestimates 
adsorption speed or equivalently slightly overestimates collector dimension. However, the 
difference would be practically imperceptible it the 9(t) dependence on dimensionless time 
was analyzed as in [1(3, [jjj]. For the sparsest collectors the relation (TlO|) is not valid and, 
therefore, parameter a cannot be determined. This is analogous to the case of Cantor Dust 
and some Cantor sets discussed in [121 ] . 



t 1 1 1 1 • r 




i i i i i i u 

2.4 2.6 2.8 3 

d 



FIG. 2: Feder's law validity. Dots represent exponents a whereas solid line is determined 
using Feder's law <Q. Inset shows power law dependence of the number of packed spheres 

on dimensionless time ( flOl) . 

The adsorption speed d6/dt may also be analyzed as a function of 9; then, it is called 
Available Surface Function (ASF). Such approach is commonly used in material science 
and especially in low coverage limit, because it allows to determine viral coefficients of the 
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equilibrium coverage state [28l. |29|. 



ASF(fl) = ^(0) = 1-^9 + C 2 9 2 + o(6 2 ), (11) 

where C\ and C% are expansion constants and they are connected with volume blocked by 
particle and a complex of two particles, respectively. For example, in a low coverage limit, 
a single ball covers v(d), but, from the point of view of another ball center, it blocks all the 
volume inside a sphere of 2r radius (see Figj3j). Therefore, using (j2J) 





FIG. 3: Left: excluded volume by single hypersphere (CI). Right: overlap of excluded 

volumes of two different particles. 

Ci(rf) = ^ = 2 d . (12) 



r 



o 



For higher coverages the exclusion volumes start to overlap, which slightly reduces the 



total amount of blocked space. Therefore, by analogy with the reasoning presented in (29J, 
parameter C% is given by 

C 2 (d) = ^ V 2 (r, d) ■ s{r, d)dr, (13) 

where 

d-l x j „™ -l _L 



V 2 {r,d) = 2 " ' v "' u; / " u sin d t dt (14) 



7i^ ± (2r ) d r CM ^ 
r[(d+l)/2]7 

is the overlapped volume and s(r, d) is (d — l)-surface of <i-sphere of radius r. Values of C\ 
and C2 measured in numerical experiments are shown in FigJU Values of C\ fits well the 
theoretical approximation. In case of C 2 the match is slightly worse, however, numerical 
data are generally arranged along the analytically obtained values. 



B. The Maximal Random Coverage Ratio 



Despite the lack of analitical results for maximal random coverage ratio for d > 1 there 
are several studies providing lower and upper limit of its value, e.g. (30(. The best relation 

8 




FIG. 4: Dependence of blocking parameters C\ and C% on collector dimension. Circles (Ci) 
and squares [Ci) represent obtained data while solid lines are analitical fits. 



describing 
proposed in 



dependence of saturated coverages on collector dimension, given by ([3]) has been 
121 ] and fits our data, see Figj5l For convenience, all numerical values 
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FIG. 5: Maximal random coverage ratio in different dimensions. Dots represents data 



obtained here and also in 
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12| and solid line is the fit ((3]). Inset shows data for 



d G [2.0; 3.0] obtained in this study. 



obtained in this study and presented on above figures are put together in TabJITl 
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G» ''max 


a C\ C2 


1.4650 0.543 ± 0.010 


2.449 0.512 


2.0 0.563 ± 0.005 


3.485 3.268 


2.4022 0.485 ± 0.005 


-0.378 ±0.012 6.105 10.194 


2.5789 0.432 ± 0.005 


-0.370 ± 0.010 6.211 9.849 


2.7268 0.415 ±0.005 


-0.348 ±0.007 7.239 14.197 


2.8540 0.398 ± 0.005 


-0.326 ±0.006 7.754 16.218 


2.9656 0.382 ± 0.005 


-0.307 ±0.004 7.945 16.899 



TABLE II: Maximal random coverages ratio, RSA kinetics and coefficient C\ and C2 were 
(fitted for 9/9 max < 0.3) obtained from RSA simulation. 

C. Pair correlation function 

Two-point correlation functions for all fractal collectors studied are presented in FigJHl 
For d close to 3, density autocorrelation function has got a typical shape. Its decay is fast 
and it has logarithmic singularity at kissing limit (r — > 2r$ ), exactly as predicted by (J4]) 
and ([5]). Situation changes in two cases in which d < 2. Then, the correlations have long 
range structure. This is due to relatively sparse collectors which strongly limit space where 
spheres could be adsorbed. The same causes unpredictable properties at small distances, 
because they do depend mainly on a specific collector structure. 



IV. SUMMARY 

Random Sequential Adsorption on Menger sponge-like fractal collectors has been stud- 
ied. For such systems, Feder's law fl6]) slightly overestimates collector dimension. Analysis 
of two-point correlation function shows that, in general, it has logarithmic singularity at 
kissing limit and its decay seems to be superexponential. The only exceptions there are the 
sparsest collectors where internal structure overcomes any universal properties of adsorption. 
Obtained maximal random coverage ratios fit well phenomenological relation proposed in 
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(e) 



(f) 



-d = 2.966 




(g) 



FIG. 6: Two-point correlation function for collectors of different dimensions. Insets show 
asymptotic character for molecules in close proximity to each other with a fit given by (jSJ). 
(a) d = 1.4650, (b) d = 2.0, (c) d = 2.4022, (d) d = 2.5789, (e) d = 2.7268 (Menger 

sponge), (f) d = 2.8540, (g) d = 2.9656 
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